Hydrodynamic theory of supersolids: Variational principle and effective Lagrangian 
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We develop an effective low-energy, long-wavelength theory of a bulk supersolid — a putative phase 
of matter with simultaneous crystallinity and Bose condensation. Using conservation laws and gen- 
eral symmetry arguments we derive an effective action that correctly describes the coupling between 
the Bose condensation and the elasticity of the solid. We use our effective action to calculate the 
correlation and response functions for the supersolid, and we show that the onset of supersolidity 
produces peaks in the response function, corresponding to propagating second sound modes in the 
solid. Throughout our work we make connections to existing work on effective theories of superflu- 
ids and normal solids, and we underscore the importance of conservation laws and symmetries in 
determining the number and character of the collective modes. 
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I. INTRODUCTION 

In 1969 Andrccv and Lifshitz proposed a novel phase of 
matter in quantum Bose crystals, wherein a Bose conden- 
sate of point defects would coexist with the crystallinity 
of the solid. 1 This is perhaps the most conceptually clear 
picture of what we now call a "supersolid," although sug- 
gestions of the coexistence (or non-coexistence) of Bose 
condensation and crystallinity can be traced to the ear- 
lier work of Penrose and Onsager 2 and Chester. 3 An- 
dreev and Lifshitz provided an elegant (albeit incom- 
plete) formulation of the hydrodynamics of supersolids, 
and predicted propagating modes analogous to second (or 
fourth) sound in liquid 4 He. Their hydrodynamic formu- 
lation was further extended by Saslow 4 and by Liu. 5 Ex- 
perimental searches for signatures of the supersolid phase 
proved fruitless 6 until recently, when Kim and Chan ob- 
served rotational inertia anomalies in solid 4 He that they 
interpreted as evidence for supersolidity. 7-9 Their work 
fueled extensive searches for further evidence of this elu- 
sive phase of matter, 10-17 and there are now a number 
of extensive reviews of the experimental and theoretical 
progress in this area — see Refs. 18, 19, 20, and 21. 

The present work is a detailed — and we believe novel — 
study of the hydrodynamics of bulk supersolids, of the 
type originally proposed by Andreev and Lifshitz. Our 
work uses conservation laws and general symmetry prin- 
ciples to derive an effective action for a supersolid. We 
rely extensively on a variational principle used to obtain 
the dynamic equations in various continuum systems: 
normal fluids, 22-24 superfluids, 24-29 normal solids, 22 ' 30 
liquid crystals, 31 trapped superfluid gases, 32 and rela- 
tivistic fluids. 33 This effective action is a powerful tool 
for calculating and elucidating the collective modes of the 
supersolid phase, and one of our important new results is 
a calculation of the correlation and response functions in 
the supersolid phase. Moving beyond linearized hydrody- 
namics, the effective action can also be use to study the 
dynamics and interaction of topological defects — vortices 
and dislocations — in the supersolid, the topics of subse- 
quent publications. 34 We should also state what this work 



is not — it is not an explanation of the recent experiments 
on possible supersolidity in 4 He, as the prevailing wisdom 
suggests that structural disorder plays a key role in most 
of the experiments, and our simplified model assumes an 
ordered solid. 

This work is organized as follows. In Sec. II we de- 
rive the supersolid hydrodynamics and the effective La- 
grangian density using the variational principle. We show 
that the equation of motion are equivalent (up to a term 
nonlinear in the elastic strain) to the hydrodynamic equa- 
tions of motion derived by Andreev and Lifshitz 1 . We 
also discuss the connection to the work of Saslow, 4 Liu, 5 
and Son. 35 In Sec. Ill we use a quadratic version of the 
Lagrangian to investigate the linearized hydrodynamics 
of a supersolid. Finally, in Sec. IV the collective modes 
and the density-density correlation function of a model 
supersolid are calculated in detail. The Appendices pro- 
vide additional detail, as an aid to the reader. 



II. VARIATIONAL PRINCIPLE AND AN 
EFFECTIVE LAGRANGIAN OF SUPERSOLIDS 

We start with a Lagrangian density for a supersolid 
in the Eulcrian description (in which all quantities are 
depicted at fixed position x and time t), 



£ss = -p sij v si v sj + ^(p6 lJ - p sij )v ni v nj 



(1) 



where p si j is the superfluid density tensor, p is the total 
density, v s is the velocity of the super-components, v n is 
the velocity of the normal-components, s is the entropy 
density, and 



RiA = '<), l> , 



(2) 



is the deformation tensor, with R the coordinate affixed 
to material elements (di = d/dxi and dt = d/dt 'm what 
follows). The first two terms in the Lagrangian density 
are the kinetic energy densities of the super-component 
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and the normal-component, respectively, and the third 
term is the internal energy density which is a function 
of p, p si j, s, and In contrast to a superfluid, Rij 

appears explicitly in Uss for a supersolid, a reflection 
of the solid's broken translational symmetry. As shown 
in Appendix B, the internal energy density satisfies the 
thermodynamic relation 



dU ss = Tds + 



I, 

P + ^{Vni 



dp - XikdRik 



-^Ki ~ V si )(v nj - V sj )dp sij , 



(3) 



where p is the chemical potential per unit mass, and 
Xij the stress tensor. Given the Lagrangian density in 
Eq. (1), the action is 



S- 



ss 



dt / d 3 x£ ss . 



(4) 



The equations of motion for a supersolid are obtained 
from variations of Sgg with respect to the dynamical 
variables. However, as illustrated in Appendix A, the 
dynamical variables are not independent and one must 
insure that conservation laws and broken symmetries are 
incorporated in the action through the use of auxiliary 
fields (Lagrange multipliers). For a three dimensional 
supersolid there are five conserved quantities: the mass, 
the entropy and the three components of the momentum. 
Among these constraints we impose only the mass and 
entropy conservation laws, and show below that the mo- 
mentum conservation is the byproduct of the variational 
principle. Conservation of mass is expressed through the 
equation of continuity, 



dtp + diji = 0, 
where the mass current ji is 

3i = Psij%j + {P$i 

The entropy conservation law is 



(5) 



(6) 



by using the Lagrange multipliers a, <f>, and with the 
result: 



1 



1 



ss 



2 PsijV Si V Bj + -(pSij - p aij )v ni v aj 



-U S s(p,Psij,s,Rij) + a 



d t s + di{sv ni ) 



+<p{ dtp + di 



Psij V sj + {PSij ~ Psij)Vnj 



dt(sRi) + dj{sRiV nj ) 



(9) 



Note that in our formulation Lin's constraint is combined 
with the entropy conservation law. 

We are now in a position to derive the hydrodynamic 
equations of motion for supcrsolids. First of all, the vari- 
ation of the action with respect to v si produces 



(10) 



Therefore, the superfluid component of the velocity is a 
potential flow, as expected [rotational flow can be ob- 
tained by introducing another constraint; see Ref. 27]. 
The remaining equations of motion are 



Sp 



dUss 
dp 



VniVs 



SPsii : 



Ss: 



dUss 1, w x 

^ = -^{Vst - v ni ){v Bj - v nj ), 



D n a D n k dU ss 

*M — ^ 1 ^ = L), 



Dt 



Dt 



ds 



(11) 



(12) 



(13) 



d t s + di(sv ai ) = 0, 



(7) 



Sv n 



in which only v n is involved because the entropy is 
transported by the normal component. Finally, we ac- 
count for the broken translational symmetry using Lin's 



constraint. 



24 



Ota + RjdiPj = -{pSij - p Hj )(v nj - v sj ), 



• SR,, 



(14) 



D n Ri 
Dt 



= 0, 



(8) 



where D n /Dt = dt + v ni di. This constraint states that 
the Lagrangian coordinates (i.e., the initial positions of 
particles) do not change along the paths of the normal 
component. Indeed, Lin's constraint was first introduced 
to generate vorticity in the Lagrangian description of an 
isentropic normal fluid. 24 We incorporate all of the con- 
straints, Eqs. (5)-(8), into the Lagrangian density Eq. (1) 



-PnA l fl ( dU SS 

Dt s j \ dRji 



0. 



(15) 



In the above equations of motion we have eliminated the 
gradient of <j> by using Eq. (10) . In addition to the derived 
equations of motion, the variations with respect to the 
Lagrange multipliers reproduce the imposed constraints, 
Eqs. (5)-(8). Therefore, Eqs. (5)-(8), (10)-(15) are the 
hydrodynamic equations for supcrsolids. 
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In the following we demonstrate that the equations 
of motion derived above are equivalent to the non- 
dissipative supersolid hydrodynamics developed by An- 
dreev and Lifshitz, 1 Saslow, 4 and Liu. 5 First, the taking 
the gradient of Eq. (11) produces the Josephson equation 



d t v s 



-difi 



1 



(16) 



where we have used the thermodynamic relation for 
dU ss /dp given in Eq. (3). Second, we derive the momen- 
tum conservation equation; the following identity simpli- 
fies the derivation: 



D(adib) Da 

-^tT = dib m 



( J^l J - adjbdtVj, (17) 



where D/Dt = d t + v t d % . Take D n /Dt of Eq. (14), and 
eliminate the Lagrange multipliers by using Eqs. (7), (8), 
(13)-(15). The result is 



Dn 

Dt 



diRjdk 



{pS i:j - p sij )(v nj 
dU ss 



= —sdi 



fdUi 



ss 



\ ds 

~ {P$ij ~ PHj)( v nj ~ V sj )d k V nk 



dR kj 

{pS jk - p sjk )(Vnk ~ V sk )diV nj 



(18) 



Third, combine Eq. (18) with the thermodynamic rela- 
tion, Eq. (3), the continuity equation, Eq. (5), and the 
Josephson equation, Eq. (16). After some algebra, we 
obtain the momentum conservation law 



dtji + djUij = 0, 



(19) 



where ji is the mass current given in Eq. (6), and LI^ is 
the (non-dissipative) stress tensor 



II,, = pv si v sj + v si pj + v nj pi - R ik Xjk 



e-Ts - fip- (v nj - v sj )pj 



(20) 



where pi = (pSij — p si j)(v n j — v S j) and e satisfies a ther- 
modynamic relation given by Eq. (B3). Note that the 
Josephson equation, Eq. (16), and the momentum con- 
servation equation, Eq. (19), are Eqs. (9) and (12) of 
Andreev and Lifshitz 1 [Andreev and Lifshitz neglected 
nonlinear strain terms, effectively replacing Rn- by Si k 
in the last term of Eq. (20) above]. Moreover, the mo- 
mentum conservation equation is equivalent to Eq. (4.16) 
of Saslow 4 when v s is taken as a Galilean velocity, and 
Eq. (3.40) of Liu 5 in the case where the superthermal 
current vanishes. 

The Lagrangian density used to derive the hydrody- 
namics of supersolids, Eq. (9), can be recast into a more 
familiar and compact form by using the equations of mo- 
tion, as illustrated for an ideal fluid in Appendix A. To 
see this, we integrate the terms involving the Lagrange 
multipliers by parts (neglecting boundary terms), and 



use Eqs. (10) and (13) to eliminate a and /3,. We then 
obtain 

^SS = ~pdt<t> ~ \p S ijdi4>dj(l) + \{p5ij ~ p Sij )VniV n j 

-{phj ~ Ps ij )v nj d l (f) - f(p, p SlJ ,T, Rij), (21) 
where / = C/gg — Ts satisfies the thermodynamic relation 



df = -sdT 
1 



dp - X lk dR. 



(Vni ~ V si )(v nj - V sj )dp s 



(22) 



When cast in this form, we see that the coupling be- 
tween the supcrfluid and the normal fluid [the fourth 
term in Eq. (21)] is — (pSij — p si j)v ni v S j; this is a "current- 
current" interaction, where the coupling constant is the 
normal fluid density. This coupling coefficient is univer- 
sal-it is determined by conservation laws and Galilean 
invariance. 

As mentioned earlier, several other authors have re- 
cently proposed Lagrangian descriptions for supersolids. 
Son 35 used symmetry-based arguments to derive an ef- 
fective Lagrangian for a supersolid. To connect to Son's 
results, we first invert Lin's constraint, Eq. (8), to obtain 



= -R^ d t Rj, 



(23) 



where R-^ = dxi/dRj and RijR-^ = Si k . Substituting 
this into Eq. (21), we obtain a Lagrangian similar in form 
to Eq. (23) of Son's paper (however, our energy density / 
depends upon p, p s , and T in addition to Rij). A differ- 
ent approach was used by Josserand et al, who applied 
a homogenization procedure to a nonlocal version of the 
Gross-Pitaevskii equation to systematically derive a long- 
wavelength Lagrangian for a supersolid. On the whole, 
our Eq. (21) agrees with their Eq. (4), once we identify 
their p(n) with our normal fluid density pS^ — p s %j- Fi- 
nally, Ye 37 proposed a phenomenological supersolid La- 
grangian; however, the Lagrangian in his Eq. (1) is not 
manifestly Galilean invariant, so that his coupling con- 
stants a Q( 3 arc arbitrary. His approach also misses certain 
nonlinear terms that are important when treating topo- 
logical defects. 



III. QUADRATIC LAGRANGIAN DENSITY 
AND THE LINEARIZED HYDRODYNAMICS OF 
SUPERSOLIDS 

In this section we discuss the propagation of collec- 
tive modes in supersolids by examining the response to 
small fluctuations away from equilibrium. The number 
of collective hydrodynamic modes of a system can be 
inferred by enumerating the system's conservation laws 
and broken symmetries. 38 Since we are more interested 
in the effect of density or defect fluctuations than ther- 
mal fluctuations, for simplicity we ignore thermal fluctu- 
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ations in what follows. For a three dimensional super- 
solid there are conservation laws for mass, three compo- 
nents of momentum, and energy; however, by ignoring 
thermal fluctuations we can omit the energy conserva- 
tion law. In addition the conservation laws, there are 
three broken translational symmetries (due to the crys- 
tallinity) and one broken gauge symmetry (due to the 
Bose- Einstein condensation). Thus, a three dimensional 
supcrsolid without thermal fluctuations has eight con- 
servation laws and broken symmetries (nine, if conser- 
vation of energy is included). Correspondingly, there 
are eight hydrodynamic modes: two pairs of the ordi- 
nary transverse propagating modes, a pair of longitudi- 
nal first sound modes, and a pair of longitudinal second 
sound modes (note that a solution of the hydrodynamic 
equations with dispersion u> = ±ck would count as two 
modes-a pair of modes, with one propagating and a sec- 
ond counter-propagating). The appearance of the longi- 
tudinal second sound modes is one of the key signatures 
of a supersolid. 

We start by establishing the notation for the small fluc- 
tuations away from equilibrium. The equilibrium value of 
the densities will be denoted with a subscript of '0', and 
the density fluctuations will be denoted by Sp, so that 
P = Pa + dp and p si j = p SQi j + Sp si y For the lattice fluc- 
tuations, let u denote the (small) deformation field away 
from the unstrained solid (i.e. the difference between the 
Lagrangian and the Eulerian coordinates): 



x = R + u. 



Then the deformation tensor becomes 



Rij = $ij - " 



%3l 



(24) 



(25) 



where Wij = diUj 



Finally, the velocity of the normal 
component is obtained by linearizing the inverted Lin's 
constraint, Eq. (23), so that to lowest order in the strain 
field the normal velocity is the time derivative of the dis- 
placement field, 



dm 



(26) 



Expanding the Lagrangian density, Eq. (21), up to second 
order in the small quantities Sp, Sp si j, 
we obtain 



did) and df 



C 



quad 

ss 



p d t 4> - Aoij-Wy - p Sp - Spd t (f> 



Po(d l (j>f 



1 
2 
1 

2 ri 
1 d\ 



dp 



dwi 



Spwij - 



ldp 
2~d~p 



Pwa - di4>) (d t Uj 

WijWlk, 



drf) 



i j 



2 dwik 



(27) 



where Pnoij = Potj — Psoij an d the thermodynamic re- 
lation for /, Eq. (22), is used. In the above expansion 
we have dropped constants which do not contribute to 



the equations of motion, and have neglected terms pro- 
portional to df /dpsij because they are of higher order 
[see Eq. (12)]; consequently, the quadratic Lagrangian 
turns out to be independent of fluctuations of the su- 
pcrfluid density. However, we have kept the first two 
terms in Eq. (27); although they are total derivatives 
and would seem to be irrelevant to the equations of mo- 
tion, they are non-trivial for topological defects such as 
vortices or dislocations. In fact, one can show 34 that the 
first term produces the Magnus force on a vortex 39 and 
the second term generates the Peach-Koehler force on 
a dislocation. 40,41 We will defer the discussion of these 
effects to a subsequent publication. 34 

Now we are in a position to study the hydrodynamic 
modes of a supersolid. The quadratic Lagrangian den- 
sity, Eq. (27), produces three linearized equations of mo- 
tions which are equivalent to the continuity equation, 
the Josephson equation and the momentum conserva- 
tion equation. Before proceeding further, it is useful to 
rewrite the quadratic Lagrangian density in terms of the 
defect density fluctuation by using one of the equations 
of motion. By varying the action with respect to 5p we 
obtain 



Sp 



dw„ 



W; 



dp 

dp 



Sp, 



(28) 



dx 


dz 


dx 


dy 


z ~ d y 





where we have used the linearized Josephson equation 
(dt4> = — Po — Sp) and the identity 



(29) 



From Eq. (28) it is clear that the density fluctuation is 
an independent hydrodynamic variable-it is not slaved 
to the lattice deformation, as would be the case for a 
commensurate solid, where dp = —po^ ■ u. Indeed, in a 
real (incommensurate) crystal a density fluctuation can 
be produced by lattice deformations or by point defects 
(vacancies and interstitials) . To highlight the role of de- 
fects we will introduce the defect density fluctuation SpA 
as our independent hydrodynamic variable, instead of 5p. 
The local defect density is defined as the difference be- 
tween the density of vacancies, py, and the density of 
interstitials, p\\ 



PA= Pi- Py- 



(30) 



The minus sign is necessary so that the total defect den- 
sity is conserved-in the bulk of the crystal vacancies and 
interstitials are created and destroyed in pairs (ignoring 
surface effects). Then we have 



Sp = 



dp 



dwi 



dp 



/'A 



dp 



A 



and from Eq. (28) we obtain 
dp 



S P : 



dwi 



dp 



dpi 



Spi 



Spi 



(31) 



(32) 
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where we have used Eq. (29) and the identity 



dx 




dx 


dx 


dz 


dy 





dy 


z " Z 


y dy 



(33) 



Equation (32) shows that a density fluctuation in an 
isothermal supersolid is caused either by a lattice defor- 
mation or by a defect density fluctuation 5p&, just as in a 
normal solid. 43-45 Following Zippelius et al. (ZHN), 44 we 

= 1. 



can identify (dp/dwij) PA 
We finally obtain 



-p Sij and (dp/dpA) 



Sp = -pawn + 5 pa, 



(34) 



which illustrates the roles of the lattice deformation wu = 
V • u and the defect density fluctuation in determining 
the total density fluctuation. We note in passing that in 
a higher order expansion of the Lagrangian density the 
terms proportional to the superfluid density fluctuation 
must also be retained in Eq. (34). This would resemble 
the "three-fluid" scenario proposed by Saslow 46 in which 
the lattice density and velocity are introduced for the 
third fluid component. 

We can now use Eq. (34) to rewrite the quadratic La- 
grangian density in terms of the defect density, with the 
result 



quad 

ss 



1 d f-i 

Pow li d t 9 - PnoijdtUid-jd - ^PQ-Qp 
-Sp A d t e - ] 2 P-.„, d,()'),0 - ^ 



Sp A Wi 



+Po 



dp 



dpA 
1 dX 



wudpi 



1 dp 



2 dp a 



ji 



2 dwi k 



dp 

WijWlk + P0g^~ 



SPA 2 

WijWkk 



(35) 



where we introduced 8 = <p + p$t. Next, we derive the 
linearized equations of motion from the Lagrangian den- 
sity. First, note that the variation with respect to 5pA 
reproduces Eq. (31) because d t is —Sp. Second, taking 
the variation with respect to 9 produces the linearized 
equation of continuity, expressed in terms of SpA ■ 



4 A 



0. 



(36) 



d t dpA + ddi 

where the defect current density is given by 

it = Psoij (djO - d t u 3 ). (37) 

We see that the defect current arises from countcrflow 
between the superfluid velocity V# and the normal fluid 
velocity dtii, and vanishes when p s oij = 0, in the normal 
state. In other words, dtSpA = in the normal state, 
in agreement with the non-dissipative description of nor- 
mal solids 44 in which defect currents are only produced 
through diffusion (i.e., the defect current is dissipative in 



the normal solid). The last equation of motion derived 
from the variation of u, is 



«2 / d P 



PnOij 



d\ 



J i 



dw 



dp 



13 dw 



ik 



dp 

dpA 

- Po 



djSpA 



dp 



dwi 



djWik = 0. 
(38) 



When the time derivative of Eq. (36) is combined with 
Eq. (31), we obtain 



dfSpA - p SQl 



dp 



dpA 



didjSpA 



p SOij didluj - p SOl 



dp 



dw 



a- 



(39) 



did.wiu = 0. 



Our linearized equations of motion, Eqs. (38) and (39), 
are equivalent to Eq. (19) of Andreev and Lifshitz. 1 

In the particular case in which the lattice sites are fixed 
(so that u = 0) we recover from Eq. (39) the fourth sound 
modes obtained by Andreev and Lifshitz, 1 which have the 
dispersion relation 



2 d P 

w =Ps 0i j 



dp 



(40) 



where we have used {d/dpA)wij = (d/dp) Wij . On the 
other hand, when there are no defect fluctuations (5pA = 
0), Eqs. (38) and (39) are combined into 



PodtUi 



dX jt 



dw 



ik 



djWik - po 



dp 



dw 



ik 



d t wi k 



-Po 



dp 



dw.; 



djWkk- 



(41) 



A mode analysis of this equation would produce six sound 
modes of an anisotropic normal solid. We see that with- 
out defects there are no additional sound modes, as ex- 
pected. 



IV. DENSITY-DENSITY CORRELATION 
FUNCTION OF A MODEL SUPERSOLID 

In this section we will calculate the density-density 
correlation function of a model supersolid, a measurable 
quantity in a light scattering experiment. However, be- 
fore delving into the calculations for a supersolid let's first 
review what's revealed by scattering light from a struc- 
tureless, normal fluid [for example, see Ref. 47]. The 
mode counting for the fluid is simple-there are five col- 
lective modes, due to conservation of mass, energy, and 
three components of momentum (in three dimensions). 
The five collective modes are a pair of transverse mo- 
mentum diffusion modes and three longitudinal modes: 



G 



a pair of propagating sound modes and a thermal diffu- 
sion mode. The density fluctuations important for light 
scattering only couple to the longitudinal modes, so three 
modes are observed: the diffusion mode appears as the 
Rayleigh peak oj = and the pair of sound modes as the 
Brillouin doublet at lu = ±cq (with a sound speed c). In 
the absence of dissipation these peaks are (5-functions; 
dissipation turns each <5-function into a Lorentzian of 
width Dq 2 , with D being an attenuation coefficient. 

What happens in a supcrfluid? In addition to the five 
conserved densities that exist in a normal fluid there is 
a broken gauge symmetry, so from mode counting we 
conclude there are six collective modes. Two of these 
are transverse momentum diffusion modes (just as for 
the normal fluid), leaving four longitudinal modes for 
the superfluid: a pair of propagating first sound modes, 
and a new pair of propagating second sound modes. In 
essence, the central Rayleigh peak in the normal fluid 
splits into a new Brillouin doublet upon passing into the 
superfluid phase. This remarkable phenomenon has been 
observed in light scattering experiments on 4 He. 48,49 We 
show below that an analogous splitting occurs in a su- 
pcrsolid, and should be observable in a light scattering 
experiment. 



fying assumptions: the solid is isotropic, and two dimen- 
sional. The isotropy causes the transverse and longitu- 
dinal modes to neatly decouple; the two dimensionality 
results in only one pair of propagating transverse modes, 
rather than two pair. Since we're interested in longi- 
tudinal fluctuations, the latter simplification is of little 
consequence to the main results of this section. With 
these assumptions, the thermodynamic relations are 



dw 



Ik 



^Sji$ik + p(SuS jk + SikSji), (42) 



dp 



7%, 



(43) 



dp 



dpi 



dp 
dp 



9 ' 

PoX 



(44) 



A. Dynamics of supersolid without dissipation 

To facilitate the calculation of the density-density cor- 
relation function for a supersolid we'll make two simpli- 



where x is the isothermal compressibility at constant 
strain, 7 is a phenomenological coupling constant be- 
tween the strain and the density, and A and p are the bare 
Lame coefficients at constant density. Then in Fourier 
space the Lagrangian density, Eq. (35), reduces to 



*pa(-Q) 

£ss = £ ( WQ) 0(Q) «l(Q) ) A I 0(-Q) 



(Pno^l + PI 2 ) Wt(Q)Wt(-Q)i 



ul(-Q 

where u) n = iu), Q = (q,w„), ul = (q • u)/g with q = |q|, Ux = u — (ui/g)q, and 



(45) 



A = 



(46) 



r 



where A = A + 2/L The collective modes are determined Setting Aa = 0, we find the longitudinal first sound 
from the determinant Aa of A: speed cl and second sound speed C2: 



Pno^n 



\PoX 



X 



1 ^ " XPo/ qi 



2 2 
1 u n 



(47) 




-2 7 

2 V \Pn POX J Pn 



xpI 7 



(48) 



7 



A 



1 



2 Pn 2 PoX 



(—- 


1 






\PnO 


PoX 


/ Pno 


.XPo 



(49) 

In particular, when p s0 = (normal solids), c 2 vanishes, 
and we only have 



4s = (A + 2/i + l/x)l Po — 27, 



(50) 



which agrees with the longitudinal sound speed obtained 
by Zippclius et al. 44 once we identify A = A ZHN +27 ZHN + 
1/X and 7 = (7 ZHN + l/x)/po- Moreover, when the Lame 
coefficients and the coupling constant 7 vanish we recover 
the sound speed of a normal fluid. As discussed earlier, 
there is one pair of transverse sound modes with speed 



c T 




(51) 



Finally, we can calculate the correlation functions from 
Eq. (45): 



(6pa(Q)Spa{-Q)) = p so q 



2 Po^ + (A-2p 7 + I/x)<z 2 



(<Va(QK(-Q)) = iq 



and 



(ul(qK(-q)> = 

PqX 



Aa 

Psq Paul - (Pol - l/x)g 2 
Po A A 

1 PoXu 2 n + Ps Q q 2 



(52) 
(53) 

(54) 



Since the density fluctuation is related to the defect den- 
sity fluctuation and the strain tensor by Eq. (34), the 
density-density correlation function becomes 



(Sp(Q)Sp(-Q)) = A 
+ B 



1 1 

iuj - c L q iuj + c L q 
1 I 



where 



A = -i 



B = —1 



iuj — c 2 q iuj + c 2 q 
PoPngcj - p SQ X 

2c L p n0 (c 2 L -c|)' 

PoPnpcl - p s0 A 
2c 2 Pn (c| -C 2 L )' 



,(55) 

(56) 
(57) 



Then, by performing the analytic continuation iu> n — 
uj + iS, the density-density response function can be ob- 
tained from the imaginary part of the density-density 
correlation function: 



Xp P (q>w) 



-ir A 
-ttB 



S(uj - c L q) - S(uj + c L q) 
5{lo - c 2 q) - 5(uj + c 2 q) 



(58) 



where we have used the identity 



= P— - — + irrS (w - u') . (59) 



uj — lo — ie uj — uj 



It is easy to show that the response function satisfies the 
thermodynamic sum rule (for the derivation of the static 
correlation function see Appendix B) 



A - ph 2 X 



and the f-sum rule 

P0Q dw 

, 7T 



ux'p P {<i,u) = poq 2 



(60) 



(61) 



B. Dynamics of supersolid with dissipation 

We continue our discussion of the density correlation 
and response functions by including dissipative terms in 
the equations of motion. As mentioned above, the dis- 
sipative terms will broaden the 5-function peaks in the 
density response function. In addition, as noted by Mar- 
tin et al. 38 , the dissipation is necessary to identify the 
"missing" defect diffusion mode in normal solids. The 
dissipative hydrodynamic equations of motion for a su- 
persolid were first obtained by Andreev and Lifshitz, 1 
who used standard entropy-production arguments to gen- 
erate the dissipative terms. For an isotropic supersolid 
(including the nonlinear term neglected by Andreev and 
Lifshitz) we have (with the new dissipative terms on the 
right hand side) 



dtp + diji = 0, 



d t ji + djUij = (did k v nk + rjd 2 v ni 



Ps(v nk - v sk ) 



(62) 



(63) 



d t Ui - v ni + v nk d k Ui + Uid k v nk = Td k X k i, (64) 



d t v si + di ( fj, + -v s 



Ps(Vnk - V s k) 



+Zdid k v nk , 



(65) 



where j = /o n v n + p s v s is the total mass current, S and A 
coefficients of viscosity, ( the bulk viscosity coefficient, -q 
the shear viscosity coefficient, T the diffusion coefficient 
for defects. 

We next linearize the dissipative hydrodynamic equa- 
tions by considering small fluctuations from the equilib- 
rium values. Writing dp and A^ in terms of Sp and dwij, 



Sp = ---Sp + jwi 
P'oX 



(66) 



5Xij = "f5 l3 Sp + XSijWkk + K w ij + w ji)- (67) 



d t u L - v n L - jTdSp - XTd 2 u L = 0, (72) 



We replace Sp and SXij into Eqs. (62)- (65), and divide 
them into transverse and longitudinal parts. The equa- 
tions of motion for transverse parts are 



Pn d t v n T ~ jld 2 u T ~ r)d 2 v n T 
where d = \[W , and 



0. 



pTd 2 ^ = 0. 



(68) 



(69) 



These equations support a propagating transverse sound 
mode with the transverse sound speed ct = \j p/pno, as 
obtained in the previous section, and an attenuation con- 
stant Ft = il + p no pT. Next, the longitudinal equations 
of motion are 



d t Sp + p SQ dv s L + p nQ dv n L = 0, 



Pn d t v n L + ^ - 7 ] 



C-2 Ps0< 7 - p s 2 A)d 2 v n L - Ps0 ad 2 v s L = 0, 



1 

pIx 



d t v B L + ^dSp+jd 2 u L -ad 2 v n L - Ps0 Ad 2 v s L = 0, (73) 



where a = S — p s0 A, and f = £ + r\. The Laplace-Fourier 
transform of Eqs. (70) - (73) yields 



(70) 



(71) where 



/ 5p(q,z) \ 

Vn L (<l,z) 

u L (q, z) 
\ v s L (q,z) J 



( Sp(q) \ 
«„ L (q) 
w L (q) 

V«s L (q),/ 



(74) 



l 



P5X 



-1 



-iz + q 2 XT 



iQPs 

Pn, 




* Pn 



9 

-9 7 



-iz + g 2 p s0 A J 



r 



(75) 



From Eq. (75) we calculate two sound speeds cl, Eq. (48), the defect diffusion mode with the diffusion constant 

D 2 = {\- p 2 xi 2 )T/p oX cl. 

We also calculate the density-density Kubo function 
rr(c£ni+n a ). from Eq. (75), 50 



and C2, Eq. (49), with two attenuation constants 

1 



Pno( c L c 2 ) 



Do 



1 



[c 2 ni + n 2 ) 



(76) 



Pn ( c L -4) 

where 

ni = ( ~ 2/? so cr + p no rA + Ps (Pn - Ps ) A , 

n 2 =-2~\ 2poPs (PoX7- l)^ + PoPn (A-PoX7 2 )r 
PoX 



(77) 



(78) 



+ PsqC + PoPsq 



Pn - Pso + Pox(A - 27p n0 ) 



A 



(79) 



Now we can see that when p s = (a normal solid), 
the second sound modes disappear but there remains 



K pp (<\,z) 



Xpp(q) iz 3 + bp p z 2 + d pp q 2 z + e pp q 2 
k B T Z 

kp>T Z ~' 



(80) 



where the static susceptibilities \pp an( i Xu L p in Eq. (80) 
are given in Appendix B, and 

Z = (z 2 - c\q 2 + izq 2 D L ){z 2 - cjq 2 + izq 2 D 2 ), (81) 

V- 

Pn 



bpp = -T 1Q 2 - PsoiPl \° Pso) Aq 2 - C 2p ^ a q 2 i (82) 



Pn 



"'pp 



X 

-i V vy, 

Pn 



-q - p s0 A7<7 + — ' q , (84) 



Pn 



Pn 
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d P u L = (A - p -/)q, 



(85) 



IPSO A 



Pn 



M 3 



+i 



PsQ 
Pn 



2vpo - C - A)A(p - 2p s0 ) 



7 g 3 .(86) 



Then the susceptibility XppCli w ) can be obtained from 
the real part of Eq. (80), 43 < 50 



Xp P (q,w) 



i<f(? L D L h{q) 



(u 2 - c 2 L q 2 ) 2 + {ujq 2 D L ) 2 



iq A c 2 D 2 I 2 (q) (lo 2 - c 2 L q 2 )I 3 (q) 

(lo 2 - c\q 2 ) 2 + (ujq 2 D 2 ) 2 (lo 2 - c\q 2 ) 2 + (uoq 2 D L ) 2 

(lo 2 - c 2 q 2 )h(q) 
(lo 2 -c 2 q 2 ) 2 + (Loq 2 D 2 ) 2 ' 

(87) 



where I\(q), -^2(9), ^3(9)1 an d -^4(9) are given in Ap- 
pendix C. 

Equation (87) is one of the central results of this paper, 
it's worth exploring some of its features and limits. First, 
one can show that Eq. (87) satisfies both the thermody- 
namic sum rule, Eq. (60), and f-sum rule, Eq. (61). Sec- 
ond, the first and second terms in Eq. (87) produce two 
Brillouin doublets centered at w = ic^g and to = ±c 2 q 
with widths D^q 2 and D 2 q 2 , respectively. The third and 
fourth terms in Eq. (87) are negligible near the Brillouin 
doublets, and in fact these terms vanish in the limit of 
zero dissipation. Therefore the non-dissipative density- 
density correlation function, Eq. (58), is obtained from 
the first two terms in the limit Dl,D 2 —> 0. Finally, for 
normal solids (p s = 0), the second term in Eq. (87) van- 
ishes, and there is only one Brillouin doublet due to the 
longitudinal first sound modes. In this case the fourth 
term in Eq. (87) becomes the Rayleigh peak of the de- 
fect diffusion mode centered at lo = 0. Therefore, we see 
that in analogy with a supcrfluid, 48 the defect diffusion 
peak that exists in a normal solid will split into a Bril- 
louin doublet of second sound modes upon entering the 
supersolid phase. 

To get a sense of the size of this effect, let's substi- 
tute some physically realistic numbers into the correla- 
tion function. Assuming p s <C po an d 7 = A = £ = 0, 
we have 



h(q) = -h(q) + 



Po 
c 2 



2i 



■ Po 

C NS 

(a -If 



Pso 



O 



a- 
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Pi 



£ 1200 - 
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-2k' 



(V C NS> 



FIG. 1: (Color online). Density-density correlation functions 
of a normal solid (black dashed line) and a supersolid (blue 
solid line). The supersolid fraction is assumed to be 10 %. 



h(q) = 



~h{q) 

(a -I) 2 po 



D A q 2 



+2q 2 



+0 



, a — 1 



n- 
Pi 



(a-l) 2 (a-3) 



a 



D A + Di 



Pso 
r 2 

(89) 



where the longitudinal sound speed of normal solid cns 
is given in Eq. (50), the defect diffusion constant -Da = 
r/x, and a = PoX c ns- We- show in Fig. 1 the normalized 
density-density correlation functions of a normal solid 
and a supersolid. We have used the first sound speed 
c ns = 550 m/s, the density p = 0.19048 g/cm 3 , the 
isothermal compressibility x = 0.29615 x 10 -8 cm s 2 /g 
for 4 He solid 51 ' 52 at the molar volume 21 cm 3 /mole, the 
viscosity of 4 Hc fluid of 2 x 10~ 5 g/cm s, the typical wave 
number involved in light scattering q^ 1 = 100 nm, and 
r = 8 x 10~ n cm 3 s/g. In Fig. 2 we show the splitting 
of the Rayleigh peak due to defect diffusion in a normal 
solid into a Brillouin doublet of second sound modes, for 
two values of the supersolid fraction. 



V. CONCLUSION 

Starting from general conservation laws and sym- 
metry principles we derived the effective action for a 
supersolid — a state of matter with simultaneous broken 
translational symmetry and Bose condensation. The re- 
sulting effective action in Eq. (21) is one of the two im- 
portant results of this work, and will be further developed 
in subsequent work on vortex and dislocation dynamics 
in supersolids. In this work, however, we used the lin- 
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FIG. 2: (Color online). Splitting of the Rayleigh peak (black 
dashed line)due to the defect diffusion mode in the normal 
solid phase into the Brillouin doublet of the second sound 
modes in the supersolid phase. The red dash-dot line is for 
ps/p = 1%, and the green solid line p s /p = 2%. 



where p is the mass density, v the velocity field, and Ujf 
the internal energy density which satisfies dUiF = \xdp. 
From the Lagrangian density we construct the action 
S = J dt J d 3 xCiF- The variational principle states that 
the equations of motion are derived from variations of the 
action with respect to all the dynamical variables. The 
naive application of this principle to the ideal fluid leads 
to the trivial equation of motion v = 0. The difficulty 
is that the dynamical variables p and v arc not indepen- 
dent, but constrained by the conservation of mass, 



d t p + d t (pvi) = 0. 



(A2) 



This constraint is incorporated into the Lagrangian den- 
sity by introducing a Lagrange multiplier <f>: 



£if = ^pv 2 - Ui F (p) + <f> 



d t p + di(pVi 



(A3) 



Then the equations of motion are obtained from varia- 
tions of the action S[p, v,</>] with respect to p, v, and 



Jp 



i 



dUiF 

= 2 V -^-Dt = °> 



D4> 



(A4) 



earized version of this action to calculate the second of 
our important results-the density-density correlation and 
response functions for a model supersolid (with isotropic 
elastic properties), see Eq. (87). In complete analogy 
with a superfiuid, we showed that the onset of superso- 
lidity causes the zero-frequency defect diffusion mode to 
split into propagating second sound modes, and from our 
calculation we can determine the spectral weight in these 
modes as well as the weight in the "normal" longitudinal 
sound waves in a solid. 
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Appendix A: Variational principle for an ideal fluid 
and the Gross-Pitaevskii action 

This appendix demonstrates the simplicity of the varia- 
tional principle for deriving the hydrodynamic equations 
of motion and the Lagrangian density for continuum sys- 
tems. Consider the simplest case of an ideal fluid (IF) 
which is irrotational, inviscid and incompressible. The 
Lagrangian density for the ideal fluid is 

Cif = \fru 2 - Ui F (p), (Al) 



— = pvi - pdi<j> = 0, (A5) 



jr = 9tP + diipVi) = 0. (A6) 

From Eq. (A5) we obtain the velocity field 

Vi = di4>, (A7) 

which implies that there is no vorticity, as expected for 
an ideal fluid. We can derive the Eulcr equation from 
Eq. (A4) by taking its gradient, using Eq. (A7), and then 
the Gibbs-Duhem relation to obtain 

P W = ~ dlP > (A8) 

where P is the pressure. The Lagrangian density may be 
cast into an equivalent form by substituting v = into 
Eq. (A3) and integrating by parts, with the result 

Cjf = ~pdt4> - ^(c^) 2 - U 1F (p). (A9) 

For comparison, consider a system of weakly interact- 
ing bosons with a condensate wave function tp(r,t) that 
satisfies Gross-Pitaevskii equation 42 

^ = -^ v V + sMV (Aio) 

This equation of motion can be derived from the La- 
grangian density 

£ GP = % w M - W] - ~ (W) ■ (VVO - 1 0A» 2 ■ 

(All) 
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Taking ip = y/ne i( P with the number density n, the Gross- 
Pitaevskii Lagrangian density becomes 

Cgp = %dtn - hnd t ^ - ^—{d in f - ^n(d t 4>) 2 - 9 -n 2 . 
2 8mn 2m 2 

(A12) 

The first term contributes ihN/2 to the action (with 
N the number of particles) and does not contribute to 
the dynamics. Identifying p = mn and <f> = (h/m)(j>, 
we see that the Gross-Pitaevskii Lagrangian density, 
Eq. (A12), is identical to the ideal fluid Lagrangian den- 
sity, Eq. (A9), with 



2m' 



•(VVP) 



2m 2 



(A13) 



Appendix B: Thermodynamic Relations and the 
Static Correlation Functions of Supersolids 



p, - v ni v si + 



^Saslow, Liu Ess and ^Saslow, Liu 

v s 2 /2. Then the differentiation of Eq. (Bl) and the use 
of Eq. (B4) give the thermodynamic relation for f/gSi 
Eq. (3). 

For a supersolid at rest we can expand the free energy 
Fgs = Egg — TS up to the second order in the density 
fluctuations bp and the strains Wij = diUj\ 

9/i r 1 d\i 



ldp, 

2~dp 



dwij 



5p Wij + - - WijWtk ■ 
(B5) 

Using Eqs. (42) - (44) for an isotropic supersolid, the free 
energy (in Fourier space) can be written as 



Fss = \m 2 u\ + - ( 5p(q) Ui(q) ) B 



where 



*p(-q) 



(B6) 



In this appendix we calculate the thermodynamic re- 
lation for the potential energy density in Eq. (1). Given 
the Lagrangian density, Eq. (1), the total energy density 
for a supersolid is defined as the sum of the kinetic energy 
densities and the internal energy density 



EsS = ~p Hj V si V sj + ^(pSij - p aij )v ui v uj 
+U S s(p, p Hj ,s,Rij). 



(Bl) 



Following Andreev and Lifshitz 1 , this total energy den- 
sity can be related to the energy density e measured in 
the frame where the super-component is at rest as 



E- 



1 



ss 



pv\ + (pSij - Psa)(Vnj ~ v aj )v ai + e, (B2) 



where e has a thermodynamic relation 



iq^ q A 



(B7) 



Then the static correlation functions can be easily read 
off from Eq. (B6): 



Xpp(q) = P (WqW-q)) 



Xu LP (q) = P (ui(q)5p(-q)) = 



- pIi 2 x) 



(B8) 



(B9) 



Appendix C: Calculation of the density-density 
correlation function 



de = Tds + pdp — X ik dR ik 

+ (l'ni ~ V si )d [(pSij - p Sij )(v nj 



(B3) 



We can obtain the thermodynamic relation for the total 
energy Ess by differentiating Eq. (B2) and using Eq. (B3) 
for de, with the result 

dE S s = Tds - \ lk dR lk - (v ni - v ai )v nj dp Hj 
1 



fi+ -(2v ni - 2v ni v si + v Bi ) 



dp 



+p aij V aj dv ai + {pSij - p SlJ )v ni dv nj . (B4) 

This thermodynamic relation agrees with Eq. (2.18) 
of Saslow 4 and Eq. (2.1) of Liu 5 after identifying 



Each term in the Kubo function given in Eq. (80) can 
be separated into the first sound part and the second 
sound part by performing a partial fraction expansion, 



g jk z 3 + b ]k z 2 + d jk q 2 z + g 2 e jk 
(z 2 — c 2 L q 2 + izD Lq 2 )(z 2 — c 2 q 2 + izD2q 2 ) 



A jk z + B 



C lk z + D 



(CI) 



z 2 — c 2 L q 2 + izD^q 2 z 2 — c\q 2 + izD 2 q 2 



where j,k = (p,u L ). Then, A jk , B jk , C jk and D jk can 
be written in terms of ajk, bjk, dj k and ej k along with 
the sound velocities (cl and c 2 ) and the attenuation co- 
efficients (Dl and D 2 ): 



A 



ik 



ajk [c 



q 2 D L (D L c 2 -D 2 c 2 L )} 



ibj k {c\D L 



D 2 c 2 L ) + d jk (c 2 L - c 2 ) + ie jk (D L - D 2 ) 



(c 



c 2 ) 2 + q 2 {D L -D 2 )(D L c 2 -D 2 c 2 L ) 



(4 



c 2 ) 2 + q 2 (D L -D 2 )(D L c 



(C2) 
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B 



ia jk c\q 2 {D L c\ ~ ggc|) + b jk c 2 L (c 2 L - ci 
(c£ - c 2 ) 2 + q 2 (D L - D 2 )(D L c 2 - D 2 c 2 



id jk cW(D L - D 2 ) + e jk [cj -4 + g 2 D L (D 2 - D L )] 

(c 2 L -c 2 ) 2 + q 2 (D L -D 2 )(D L c 2 -D 2 c 2 L ) ■ 1 > 



The coefficients Cj k and Dj k are the same as Aj k and ^(q) = Xppi^Bpp + Xu L p{<i)B pUL — iq 2 D , (C6) 

Bj k , respectively, but two sound velocities and two at- 
tenuation coefficients must be interchanged. Then the 
functions defined in the density-density correlation func- 
tion, Eq. (87), are given by 

h(<l) = X PP (q)A pp + Xu LP {q)A pUL , (C4) J 4 ( q ) = Xpp (q)D pp + X u LP (<l)Dpu L - iq 2 D 2 I 2 {q). (C7) 



h(q) = Xp P (q)C PP + Xu LP {<£)Cp UL , 



(C5) 
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